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freespace  and  the  time  delay  is  therefore  longer. 
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imple  empirical  expression  for  range  error  for  elevation  angles  above  5°.  The 
!RPL  Reference  Atmosphere  1958  was  the  model  used.  A  simple  regression 
line  that  yields  range  error  corrections  within  1  percent  of  the  values  obtained 
by  means  of  more  sophisticated  techniques  was  derived.  / 
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Satellite  systems  like  the  621B  have  been  proposed  for  air  traffic  control  and 
navigation.  Such  systems  are  designed  to  provide  aircraft  position  and  velocity 
data  based  on  time  deiay  measurements  of  propagation  between  the  aircraft  and  a 
network  of  satellites.  Since  the  index  of  refraction  of  the  troposphere  is  greater 
than  unity,  radiowaves  propagate  through  the  troposphere  slower  than  through 
freespace  and  the  time  delay  is  therefore  longer. 

Range  error  corrections  calculated  from  various  models  of  index  of  refrac¬ 
tion  generally  require  the  use  of  numerical  methods  and  a  large-capacity  high¬ 
speed  computer,  particularly  for  low  elevation  angles.  This  paper  presents  a 
simple  empirical  expression  for  range  error  for  elevation  angles  above  5°.  The 
CRPL  Reference  Atmosphere  1958  was  the  model  used.  A  simple  regression 
line  that  yields  range  error  corrections  within  1  percent  of  the  values  obtained 
by  means  of  more  sophisticated  techniques  was  derived. 
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Corrections  for  Tropospheric  Range  Error 


1.  INTRODUCTION 

Satellite  systems  like  the  621B  have  been  proposed  for  air  traffic  control  and 
navigation.  Such  systems  are  designed  to  provide  aircraft  position  and  velocity 
data  based  on  time  delay  measurements  of  propagation  between  aircraft  and  a 
network  of  satellites. 

Since  the  index  of  refraction  of  the  troposphere  is  greater  than  unity,  propa¬ 
gation  of  radiowaves  through  the  troposphere  is  slower  than  through  free  space. 
Also,  the  variation  in  the  index  of  refraction  as  a  function  of  altitude  produces  a 
bending  of  the  propagated  wave,  with  the  net  result  an  increase  in  the  path  length. 

To  determine  the  tropospheric  range  error  that  produces  the  time  delay,  the 
index  of  refraction  of  the  troposphere  along  the  ray  path  must  be  known.  A  num¬ 
ber  of  models  of  the  index  of  refraction  as  a  function  of  altitude  have  been  used 
for  calculating  the  tropospheric  range  error.  The  expressions  obtained  are 
rather  complicated,  however,  particularly  for  low  elevation  angles.  Numerical 
methods  and  a  large -capacity  high-speed  computer  must  be  resorted  to  for  a 
solution,  which  may  nevertheless  be  difficult  to  achieve  in  real  time  under  opera¬ 
tional  conditions. 

The  atmospheric  index  of  refraction  decreases  approximately  exponentially 
with  altitude,  and  so  range  error  decreases  with  increasing  elevation  angle.  For 
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angles  above  15°  the  range  error  correction  can  be  obtained  from  a  very  simple 
expression.  This  paper  presents  the  derivation  of  a  simple  but  accurate  empiri¬ 
cal  expression  for  range  error  correction  for  elevation  angles  down  to  5°. 


2.  REFRACTIVITY  MODELS 

Numerous  refractivity  (N)  profiles  have  been  produced  over  the  years,  from 
samples  measured  directly  with  a  refractometer  or  calculated  from  measure¬ 
ments  of  temperature,  pressure,  and  relative  humidity.  A  number  of  models 
have  been  generated  from  these  data  for  the  purpose  of  correcting  for  tropo¬ 
spheric  effects  on  propagation.  Historically,  the  earth*  s  radius  model  formu¬ 
lated  by  Schelleng  et  al  (1933)  was  first  used  for  line-of-sight  communications 
problems.  It  was  shown  that  by  assuming  an  earth  having  a  radius  of  about  4/3 
that  of  the  actual  earth,  radio  rays  could  be  drawn  as  straight  lines. 

The  4/3  earth  model  works  well  for  propagation  paths  at  low  altitudes,  where 
ray  paths  are  within  about  2  km  of  the  earth*  s  surface,  but  not  for  those  at  higher 
altitudes.  From  an  examination  of  many  years  of  N -profile  data  for  various  cli¬ 
mates  it  appears  that  N  decreases  approximately  exponentially  above  an  altitude 
of  about  1  km  and  that  the  variation  in  refractivity  is  minimum  at  an  altitude  of 
9  km,  with  a  range  of  about  only  8  N  units  and  an  average  value  of  104.8  N  units 
(Bean  and  Dutton,  1966,  p.  61).  The  Rocket  Panel  data  (1952),  the  ARDC  Model 
Atmosphere  1956  (1957)  and  work  by  Dubin  (1954)  show  that  above  9  km,  refrac¬ 
tivity  decreases  exponentially.  These  data  were  used  in  formulating  the  model 
known  as  the  CRPL  Reference  Atmosphere  1958,  which  includes  the  following 
expressions  for  refractivity  for  three  ranges  of  altitude: 

N(h)  as  N  +  (h  -h  )  AN,  h  <h£h+lkm, 

s  s  s  s 

where 

AN  =  -7.  32  exp(0. 005577)  Ng; 

N(h)  =  Nj  exp  [  -c(h  -hs  -  1)],  hg+  j  A.  h  9  km  , 
where 

i  .  Ni 

c  =  8^rlnTc5: 

s 

N(h)  =  105  exp  [  "0. 1424(h  -  9)],  h>9km, 
where 

N(h)  »  refractivity  at  altitude  h, 
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Nfi  =  surface  refractivity, 

Nj  =  refractivity  at  altitude  1  km  above  the  surface, 

hfi  =  height  above  sea  level. 

To  simplify  these  expressions,  the  following  single  exponential  (referred  to 
as  the  CRPL  Exponential  Reference  Atmosphere  1958)  was  generated: 

N(h)  =  Ng  exp  [-ce(h  -  hg)  ] , 
where 

Ns 

Ce  =  111  Ng+  AN  ' 

Although  this  model  is  in  good  agreement  with  refractivity  data  below  8  km,  the 
values  it  gives  at  higher  altitudes  are  too  low.  For  example,  at  an  altitude  of 
9  km  the  average  refractivity  is  about  105,  with  minimum  and  maximum  values  of 
100  and  108  respectively;  but  at  this  altitude  the  Exponential  Reference  Atmos¬ 
phere  gives  a  value  of  only  85  for  Ng  =  31  3  (average  surface  refractivity  in  the 
United  States).  This  model  is  therefore  not  expected  to  be  very  useful  for  range 
error  calculations. 

Hopfield  (1969)  noted  that  if  the  refractivity  as  a  function  of  height  is  repre¬ 
sented  by  an  exponential,  it  is  not  integrable  in  closed  form,  whereas  if  it  has 
the  form 


N  =  k(h  -hf,  h  <  h  . 
o  o 

it  is  integrable-  She  showed  that  representing  the  dry  and  wet  terms  of  refractiv¬ 
ity  by  quartic  equations  (p  =  4)  gave  good  agreement  with  range  error  and  doppler 
data  (Hopfield,  1971). 

A  number  of  other  models  of  refractivity  have  been  proposed.  For  the  most 
part  they  are  variations  of  those  mentioned  above  and  will  not  be  discussed. 

Probably  the  best  source  of  refractivity  data  is  the  CRPL  Standard  Atmos¬ 
pheric  Radio  Refractive  Index  Sample  (Bean.  Cahoon,  and  Thayer,  1960).  It 
consists  of  a  sample  of  77 N  profiles  selected  from  thirteen  radiosonde  stations 
representative  of  the  major  geographic  and  climatic  types  of  the  world.  Range 
errors  were  calculated  for  each  profile  as  a  function  of  elevation  angle,  using  ray- 

tracing;  and  least -square  regression  lines  of  the  form  AR  =  a 4*  bN  were  ob- 

o  e  s 

tained  for  some  angles.  For  elevation  angles  above  5°,  range  errors  calculated 
from  the  CRPL  Reference  Atmosphere  1958  were  found  to  be  in  good  agreement 
with  those  obtained  from  this  regression  line  (Norton,  1964).  The  CRPL  Reference 
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Atmosphere  1958  was  therefore  selected  for  the  range  error  calculations.  Fortu¬ 
nately,  results  based  on  this  model  were  available  from  Mill  man  (1970),  who  had 
generated  extensive  curves  of  range  error  as  a  function  of  elevation  angle,  surface 
refractivity,  and  height.  By  selecting  a  model  consisting  of  stratified  layers  only 
50  m  thick,  from  ground  level  to  an  altitude  of  30km,  and  taking  into  accound  re¬ 
fractive  bending,  he  was  able  to  compute  very  accurate  range  errors. 


3.  FORMULATION  OF  SIM  I’Ll  FIFO  EXPRESSION  FOR  RANGE  ERROR 


The  object  in  the  work  reported  here  was  to  generate  a  simple  expression 
that  would  provide  range  errors  in  close  agreement  with  those  computed  by 
Millman  (1970)  for  the  CRPL  Reference  Atmosphere  1958  for  elevation  angles 
above  5°.  Bean  and  Dutton  (1966,  pp.  337-339),  using  ray-tracing  techniques, 
had  shown  that  range  error  due  to  bending  of  the  propagated  wave  is  small  for 
elevation  angles  above  3°,  and  had  derived  the  following  general  expression  for 
range  error: 


They  estimated  that  the  summation  terms  contribute  about  3  percent  or  less  for 
10°,  and  about  12  percent  for  s  3°,  and  obtained  a  regression  curve  for 

10'3/  N(h) dh  =  1 . 4588  +  0.002961  IN  ,  (2) 

«/o  s 

which  is  the  limiting  form  of  the  total  range  error  based  on  the  CRPL  Standard 
Sample  for  0  =  90°. 

Combining  the  first  term  on  the  righthand  side  of  Eq.  (1)  with  Eq.  (2)  gives: 


AR  =  esc  0  (1.4588  +  0.0029611  N  )  ,  (3) 

e  o  s 

where 

AR^  =  range  error  in  meters, 

9  a  elevation  angle, 

Ns  3  surface  refractivity. 

A  comparison  of  range  errors  obtained  from  Eq.  (3)  with  those  computed  by 


Millman  (1970)  (see  Figure  1)  shows  that  the  differences  increase  approximately 
exponentially  with  decreasing  elevation  angle  for  constant  surface  refractivity. 


ELEVATION  ANGLE  (DEGREES) 

Figure  1.  Differences  Between  Approximate  and  Accurate 
Expressions  for  Range  Error  as  a  Function  of  Angle  (N 
constant)  s 

Plotted  as  a  function  of  refractivity  for  constant  elevation  angle  (Figure  2)  the  dif 
ferences  are  shown  to  be  quadratic  in  form.  Hence,  adding  a  term  of  the  form 

A(0,Ns)  =  A  [c/>.  c3(Ns-c4)2  +cj 
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to  Eq.  ( 3)  should  decrease  the  differences.  As  shown  in  Appendix  A, 

A(0,Ns)  =  jo. 00586  <Ns  -  360)2  +  294j  0-2'  30  .  (4) 

Therefore, 

AR  a  U9t  0.00072  Ife  _  r0#00586  (N  _  360)2+  294]  «’2,  30  .  (5) 

eo  sin  v  L  s  j 


where 


AR  =  range  error  in  feet, 
eo 

6  =  elevation  angle  in  degrees, 

Nfl  =  refractivity  in  N units  at  surface. 

Equation  (5)  yields  range  errors  from  the  earth' s  surface  through  the  total  tropo¬ 
sphere  that  are  within  ±0. 5  ft  of  those  computed  by  Millman  (1970). 


Figure  2.  Differences  Between  Approximate  and  Accurate 
Expressions  for  Range  Error  as  a  Function  of  Surface 
Refractivity  (0  constant) 
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The  next  step  is  to  obtain  an  expression  for  range  error  for  an  aircraft  at 
altitude  h  through  the  total  atmosphere. 

Range  error  for  a  path  from  sea  level  to  an  aircraft  at  altitude  h  has  been 
computed  by  Millman  (1961f  pp.  51-52),  from  whose  results  it  should  be  possible 
to  estimate  the  range  error  for  the  path  between  an  aircraft  and  a  satellite.  For 
example,  if  the  total  range  error  from  sea  level  to  a  satellite  is  AR ,  and  the 
range  error  from  sea  level  to  an  altitude  h  is  ARe,  then  the  range  error  AR^  for 
an  aircraft  at  altitude  h  to  the  satellite  is  approximately  AReh  =  AReQ  *  A^e' 

From  Figure  3,  where  Millman's  (1970)  range  errors  are  plotted  as  a  function  of 
altitude  for  several  elevation  angles  and  surface  refractivities  of  N  =  280  and 
N  =  360,  it  can  be  seen  that  the  height  dependence  of  range  error  is  almost  expo¬ 
nential  at  lo#  altitudes  but  falls  off  more  rapidly  at  higher  altitudes. 

In  Appendix  B  it  is  shown  that  the  range  error  from  an  aircraft  at  altitude  h 
through  the  total  troposphere  to  a  satellite  can  be  approximated  by 

AReh  =  AR  exp  -  [(6.07  X  10'5NS  +  0.021  3) h  +  -  1.58  X  10_4)h2].  (6) 

Equation  (6)  yields  range  errors  that  agree  with  Millman1  s  values  to  within  ±0.  7  ft. 
It  has  now  been  shown  that  range  errors  in  excellent  agreement  with  those  obtained 
based  on  the  CRPL  Reference  Atmosphere  1958  can  be  computed  quite  easily  from 
Eq.  (6)  for  elevation  angles  greater  than  5°.  These  have  been  tabulated  in 
Appendix  C. 

It  has  been  assumed  that  the  surface  refractivity  is  known.  The  accepted 
expression  (Bean  and  Dutton,  1966,  p.  7)  for  the  calculation  of  refractivity  from 
temperature,  pressure,  and  water  vapor  pressure  is 

N(h)  =  77.6SH+  3.73Xl0'5-^4  ,  (7) 

IX  h>  T(hr 

where 

P(h)  =  pressure  in  millibars  at  altitude  h, 

T(h)  =  temperature  in  degrees  Kelvin  at  altitude  h, 
e(h)  =  water  vapor  pressure  in  millibars  at  altitude  h. 

For  an  aircraft  at  altitude  h,  the  refractivity  is  related  to  the  surface  refractiv¬ 
ity  by 

N(h)  =  N  exp(-h/H)  , 
s 

where  H  is  ?i  scale  height  of  about  7  km,  (Bean  and  Dutton,  1966,  p.  15).  There¬ 
fore  Ng  =  N(h)  exp(0.043h),  where  h  is  in  thousands  of  feet. 
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If  the  temperature  and  pressure  are  known,  but  the  water  vapor  pressure  is 
unknown,  it  becomes  necessary  to  estimate  the  refractivity  from  average  values 
of  relative  humidity.  Bean  and  Dutton  (1966,  pp.  398-403)  have  tabulated  refrac¬ 
tivity  as  a  function  of  temperature  and  relative  humidity  for  pressures  of  1000, 
850,  and  700  millibars,  which  correspond  approximately  to  altitudes  of  0,  5000, 
and  10000  ft,  respectively.  From  the  U.S.  Standard  Atmosphere  Supplements 
(1966,  p,  12)  it  can  be  seen  that  for  all  northern  latitudes — and  presumably  for 
all  southern  latitudes— the  average  relative  humidities  at  0,  5000,  and  10000  ft 
are  respectively  78,  70,  and  50  percent. 

In  Appendix  D  the  wet  term  contribution  to  the  refractivity  is  listed  as  a 
function  of  temperature  and  pressure.  It  is  shown  that  it  can  be  approximated  by 
the  function 

Nw(h)  =  exp  jo.  0824[T(h)-232]  - 

-1.88Xl0'4[T(h)-232]2-4.5h2[T(h)-203]Xl0'5(  .  (8) 

Total  refractivity  for  the  case  where  the  relative  humidity  is  unknown  is  therefore 

N(h)  =  Nd(h)  +  Nw(h) 

=  77,6  TUi)  +  expj0.0824[T(h)-  232]- 

-1.88X10'4[T(h)-232)2-4.5h2[T(h)-203]Xl0'5|  .  (9) 

If  real-time  meteorologic  data  are  not  available,  there  are  two  choices.  The  first 
is  to  use  an  average  seasonal  refractivity  (or  temperature,  pressure,  relative 
humidity)  for  the  location  of  the  aircraft;  the  second  is  to  use  the  refractivities  in 
Table  1,  which  are  calculated  from  the  average  summer  and  winter  temperature, 
pressure,  and  relative  humidity  as  a  function  of  latitude. 


Table  1,  Average  Surface  Refractivity 


Latitude 

(degrees) 

January 

July 

0  to  22.5 

369 

369 

22.5  to  37.5 

328 

382 

37.5  to  52.5 

307 

344 

52.5  to  67.5 

310 

324 

67.5  to  82.5 

316 

315 

10 


4.  DISCUSSION  OF  ERRORS 

It  is  extremely  difficult  to  determine  the  accuracy  with  which  tropospheric 
range  error  can  be  estimated  since  it  is  not  practicable  to  measure  range  error 
directly.  Although  range  error  can  be  accurately  calculated  by  means  of  ray¬ 
tracing  techniques  if  the  refractivity  along  the  path  is  known,  it  is  not  always 
possible  to  ascertain  atmospheric  refractivity  since  it  is  variously  affected  by 
winds,  clouds,  and  precipitation.  In  addition,  the  assumption  that  the  atmosph ere 
is  horizontally  stratified  is  not  always  valid.  It  is  particularly  difficult  to  acc  u¬ 
rately  measure  refractivity  at  the  higher  altitudes.  There  is  no  uniformity  of 
sensors  used  in  radiosondes;  although  temperature  and  pressure  sensors  generally 
show  fair  agreement,  humidity  sensors  are  erratic  (Bean  and  Dutton,  1966, 
pp.  23-47).  Refractometers  are  more  capable  of  accurately  measuring  refrac- 
tivities,  but  they  are  relatively  expensive  and  complex  to  operate.  With  these 
limitations  in  mind,  range  errors  AR£  that  had  been  calculated  for  the  CRPL 
Standard  Sample  were  examined  by  Norton  (1964),  who  plotted  the  standard  devia¬ 
tion  of  ARe/Re  as  a  function  of  elevation  angle  for  ranges  of  Rg  =  10,  100,  and 
1000  km  for  cases  where  the  surface  refractivity  is  both  known  and  unknown.  The 
curves  corresponding  to  Rg  =  1000  km  were  selected  for  the  work  of  determining 
the  range  error  to  a  satellite,  and  the  corresponding  standard  deviations  of  ARg 
plotted  in  Figure  4. 

If  it  is  assumed  that  the  range  errors  obtained  by  using  the  CRPL  Reference 
Atmosphere  1958  have  approximately  the  same  standard  deviation  as  those  ob¬ 
tained  by  using  the  CRPL  Standard  Sample,  then  the  curves  in  Figure  4  can  be 
used  to  estimate  the  accuracy  of  the  range  error  corrections.  With  Ng  known, 
the  standard  deviations  range  from  1 .  3  ft  at  an  elevation  angle  of  5°  to  about 

0. 1  ft  at  90°.  With  N  unknown,  the  standard  deviations  increase  to  6.  7ft  and 
s 

0.  7  ft  respectively.  If  the  temperature  and  pressure  are  known  but  the  humidity 
is  unknown,  the  standard  deviation  is  presumed  to  lie  somewhere  between  these 
two  ranges. 

Equation  (6)  has  been  shown  to  produce  range  error  corrections  that  differ  by 
less  than  1  ft  from  the  more  accurate  values  of  Millman  (1970).  This  error  must 
also  be  taken  into  account. 

At  best,  only  a  qualitative  estimate  of  the  accuracy  of  range  error  corrections 
can  be  presented.  The  largest  uncertainty  exists  in  cases  of  low  elevation  angles. 
If  Ng  is  known,  typical  range  error  corrections  accurate  to  2  or  3ft  can  be  expect¬ 
ed;  in  rare  cases  (3a)  they  may  approach  10  ft.  If  Ng  is  unknown,  the  accuracy  of 
the  range  error  correction  decreases  significantly,  with  typical  errors  on  the 
order  of  8  ft  and  atypical  errors  approaching  20  ft.  These  accuracies  are  consis¬ 
tent  with  those  estimated  by  Freeman  (1964)  and  Evans  (1971). 
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Figure  4,  Standard  Deviation  of  Range  Error  as  a  Function  of 
Elevation  Angle 


5.  COMMENTS 

Tropospheric  refraction  produces  range  errors  approaching  100  ft  for 
ground -to -satellite  paths  at  elevation  angles  of  5°;  if  these  errors  go  uncorrected 
they  can  seriously  degrade  air  traffic  control  and  navigation  systems  based  on 
time  delay  measurements. 

For  calculating  range  error  corrections,  the  CRPL  Reference  Atmosphere 
1958  was  selected  as  the  most  suitable  of  the  several  models  of  refractivity  pro¬ 
files  that  were  reviewed.  The  complexity  of  this  model  generally  requires 
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numerical  methods  and  a  large -capacity  digital  computer  for  accurately  calculat¬ 
ing  range  error,  particularly  for  low  elevation  angles.  The  simple  regression 
lines  that  have  been  derived  yield  range  error  corrections  within  about  1  percent 
of  these  more  accurate  values  for  elevation  angles  5°  and  above,  allowing  real¬ 
time  computation  of  range  error  with  a  small  computer. 

The  true  accuracy  of  the  range  error  correction  remains  unknown,  for  reasons 
stated  in  Sec.  4.  On  the  basis  of  range  error  statistics,  however,  it  is  estimated 
that  when  the  surface  refractivity  is  known,  the  calculated  values  are  probably  close 
to  the  true  values  90  percent  of  the  time;  and  that  when  the  surface  refractivity  is 
unknown,  then  errors  on  the  order  of  10  percent  of  the  true  value  can  be  expected. 
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Appendix  A 

Correction  of  Range  Error  From  Sea  Level 

As  seen  in  Figure  1,  the  range  error  difference  A  has  an  elevation  angle  de¬ 
pendence  that  is  approximately  linear  on  a  log-log  scale.  Therefore,  A  should 
have  the  form: 

lnA(0,Ng  constant)  =  lnCj  +  c^lnO  , 
or, 

c2 

A(0fN  constant)  =  c,0 

S  A 

With  N  =  240, 

for  0  =  6°:  In  5. 9  =  lnc^  +  c^  In  6°  ; 

for  0  =  10°:  In  1 . 8  =  In  Cj  +  c^  In  10°  . 

Solving  for  Cj  and  c^  yields 

c}  =  378  ; 
c2  =  -2. 30  . 
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Therefore, 

A(0,N  =  240)  =  3780-2’30. 

s 

In  Figure  2  the  range  error  difference  was  seen  to  have  a  refractivity  de¬ 
pendence  that  is  slightly  quadratic,  with  a  minimum  at  about  360.  Therefore, 
A(Ng,  6  constant)  should  have  the  form 

AtN  ,6  constant)  =  c7(360-N  )2  +  c.  . 
s  3  s  4 

With  6  =  6°, 

for  Ng  =  280:  5.  3  =  c3(  360-280)2  +  c4  ; 

for  N  =  360:  4.  7  =  cA  . 

s  4 

Therefore, 

c3  =  9. 38  X  lo'5  , 

A(Ns,  0=6^  =  9.  38  X  10-5  (Ng  -  360)2  +  4.  7  . 

To  combine  these  results  into  a  single  expression,  noting  that  A  has  a  stronger 
^-dependence  than  N -dependence,  let  it  be  required  that 

c5[A(Ns,  0=6°)1  =  378,  for  N  =  240, 

c5[9.  38 X  10~5  (240-360)2  +  4.  7)  =  378, 

c5  a  62. 5  . 


Then 


A(0,NS)  =  [0.00586  (Ng  -360)2  4  2941  0'2  3°  . 


[(4)] 
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Appendix  B 

Correction  of  Range  Error  From  an  Aircraft 


Range  errors  plotted  versus  altitude  on  semilog  paper  (Figure  3)  are  seen  to 
be  almost  linear,  which  indicates  an  exponential  altitude  dependence.  To  com¬ 
pensate  for  the  nonlinearity  at  the  higher  altitudes,  an  additional  term  is  added. 
The  expression  for  range  error  as  a  function  of  height  is  then 

lnAReh  =  lnAReo+  c,h  +  c£Z  . 

The  term  c.  is  determined  from  values  of  AR  u  at  h  =  0  and  h  =  2  km;  c0  is 

1  #  eh  2 

determined  from  values  of  AR  .  at  h  =  0  and  h  =  20  km.  Thus, 

eh 

In  AR  .  -  lnAR_ 

_  e2  eo 

C1  =  2 

In  AR  “  In  AR  -  20  k, 

e20  eo  1 

c2  =  400 

These  have  been  computed  for  N  =  280  and  N  =  360  at  five  elevation  angles,  as 

s  s 

follows; 


Ns=280 

Ns= 

=360 

C1  c2 

C1 

C2 

6  -0.130  -0.0012 

-0.143 

-0.0007 

10  -0.126  -0.0012 

-0.143 

-0.0006 

15  -0.123  -0.0014 

-0.141 

-0. 0006 

20  -0.122  -0.0014 

-0.140 

-0.0004 

30  -0.121  -0.0012 

-0.139 

-0.0009 

On  the  basis  of  the  above  results,  the  constants 

for  Ng  =  280  are  taken  as:  =  -0. 126,  c2  =  -0.0013  ; 
for  Ng  =  360  are  taken  as:  =  -0. 142,  c2  =  -0.0006  . 

Since  these  constants  vary  slightly  with  Ng,  a  linear  dependence  is  assumed, 
as  follows: 


280  a,  +  a  =  -0. 126 
i  o 


360  a,  +  a  =-0.142 
1  o 


3lq  -  -0.07,  a^  s  -0.0002 
257+  bQ  = -0.0013 


150  +  b0  =  -0.0006 


bQ  =  0.0017,  bjSs-O.t'3 


=  -0.07  -0.0002 N 


c,  a  -  °‘83  +  0.0017 

S 


AReh  =  AReo  exp  -|(0.0002NS  +  0.07)h  +  -0.00n)h2J  , 

where  h  is  in  kilometers.  Converting  h  into  thousands  of  feet  yields 

AReh  =  AReo  exp  -j(6.07  X  10"5Ns  +  0.0213)h  +  -  1.58  X  10-^h2j.  [(6)] 
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Appendix  C 

Table  of  Range  Errors 


RANGE  ERRORS  (feet) 
for  h  •  sea  level 


20 


RANGE  ERRORS  (feet) 
for  h  -  2000  ft 


RANGE  ERRORS  (feet) 
for  h  ■  4000  ft 
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RANGE  ERRORS  (feet) 
for  h  ■  6000  ft 


RANGE  ERRORS  (feet) 
for  h  •  8000  ft 


RANGE  ERRORS  (feet) 
for  h  •  10000  ft 


RANGE  ERRORS  (feet) 
for  h  ■  15000  ft 


RANGE  ERRORS  (feet) 
for  h  -20000  ft 


27 


RANGE  ERRORS  (feet) 
for  h  ■  25000  ft 
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RANGE  ERRORS  (feet! 
for  h  -30000  ft 


RANGE  ERRORS  (feet) 
for  h  -60000  ft 


RANGE  ERRORS  (feet) 
for  h  •  80000  ft 


33 


UJ  j-  J  lO  Kt  CN  ICJIT.  OONNlAaomHOtCO  CO 

at  U)  ^  M  H  O  CONClAlA^I^NNNHHH 


W  H  H  H  »-lf“l»-lOOOOOOOOOOOOO 


id  in  ki  u>  w  cm  iD-jmooKi»niDooK»  h  cn  «o  ao 

Kl  C>  lO  J  «M  r-*  O  CONU>lAlA-?l^CNN  N  H  H  H 


04  r-*  r-4  r-4  r-4  r~*  r-i  OOOOQO  O  CD  O  O  O  O  CD 


MNKtNiANHiAnmaoNNOaomooteoeo 

K1  C»  O  J  04  r-4  O  CONUtAtA^^NNNHHH 


04  *— 4  «— 4  *— 4  «— 1  H  H  OOOOOOOOO  oooo 


OO  CO  O  C*  04  o  cn  mNKINHHlANldOOCON 

N  CO  lO  m  04  H  at  CONIOUtUt^inNNNHHH 


C4HHHHHOOOOOOOOOQOOOO 


<Jt  H  J  J  CO  lO  lO  HOHinoO^lDNOCONO 

H  CO  IA  m  H  O  at  CONUJlAUtjKtNNNHHH 


N  >H  H  H  H  H  O  OOOOOOOOOOOOO 


oi  o  id  cn  in  jiamniaci^ioo  ^  o  n  io  in  ia 

atlOKtHOatCOrslOlAJJinKIMWHHHH 


r-l«'-l*-l<-lf-IOOOOOOOOOOOOOOO 


r 
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Appendix.D 

Wet  Term  Contribution  to  Refractivity  (Relative  Humidity  Unknown) 


The  wet  term  contribution  to  refractivity  for  average  relative  humidities  of 

78,  70,  and  50  percent,  at  pressures  of  1000,  850,  and  700  mb,  respectively,  is 

obtained  by  subtracting  the  dry  term  refractivity  (zero  percent  relative  humidity) 

from  the  total  refractivity  (Bean  and  Dutton,  1966,  pp.  398-403). 

When  the  wet  term  of  the  surface  refractivity  (N  )  is  plotted  against  surface 

w  s 

temperature  on  semilog  paper  (Figure  Dl)  it  appears  that  it  can  be  approximated 
with  a  regression  line  of  the  form 

^ws^l  <T’co>+  c2(T_co)2  * 

Setting 

T  =  243°k,  In  2.  5  =  Cj  (243  -  cq)  ; 

T  =  253°K,  In  5.  7  =  Cj  (253 -cq)  ; 

we  get 

c  =  232,  c1  =  0.0824  . 
o  1 

Then  for 

T  =  303°K,  In  1  34.  6  =  0.0824(303-232)  +  c2(303-232)2  , 
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we  get 

c2  =  -1.88X  10"4  . 

At  sea  level,  therefore, 

=  exp  [0.0824<T-232)  -1.88X  10~4(T-2322]  . 

In  Figure  D2  the  wet  component  of  refractivity  is  plotted  as  a  function  of  pres¬ 
sure  for  a  series  of  constant  temperatures.  It  can  be  seen  that  this  family  of 
curves  has  a  pressure  (or  height)  dependence  of  the  form 

’U’N.s'*"*-'*2'- 

where  pressures  of  1000,  850,  and  700  mb  are  assumed  to  correspond  to  altitudes 
h  of  approximately  0,  5000,  and  10000  ft,  respectively,  and  k  is  a  constant.  It 
can  be  shown  that 

k  =  -4.5(T-203)X10'5  . 

Thus, 

=  exp  [0.0824(T-232)  -1.88X  10'4(T-232)2-4.5  h2(T-203)X  10-5]  . 


